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Abstract 

For a general Hamiltonian appropriate to a single canonical degree of freedom, we char- 
acterize and define a universal propagator with the property that it correctly evolves the 
coherent-state Hilbert space representatives for an arbitrary fiducial vector. The universal 
propagator is explicitly constructed for the harmonic oscillator, with a result that differs 
from the conventional propagators for this system. 


1 Introduction 

Canonical coherent states, and the coherent state propagator they engender, have been around 
for over three decades .* -3 In essence, their construction is simplicity itself. Let P and Q denote 
an irreducible pair of self-adjoint Heisenberg operators satisfying [Q,P] = i(h = 1)> and let 

= e~‘ qF e i ’%) 

denote a family of normalized states defined for a fixed fiducial vector | q) , (q\ q) = 1 , and for all 
(p,<?) € R 2 . These states are the canonical coherent states and they admit a resolution of unity 
in the form 

] \p,q;ri)(p,<i',i\<ipdq/ 2 * = i , 

for any 1 77 ), when integrated over all phase space . 2 These states lead to a representation of Hilbert 
space H by bounded, continuous functions, 

,(p,?) = (p,q-,v |V>), 

defined for all \tp) € H, that evidently depend on the choice of \q), although that dependence is 
often left implicit. An inner product in this representation is afforded by 

(<f>\ip) = J <t>‘ n (p, q^ip, q)dpdq /2 tt , 

an integral which removes all trace of the fiducial vector \q). 
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1.1 Propagators 

The abstract Schrodinger equation 

= W), 

involving the self-adjoint Hamiltonian H, is formally solved with the aid of the evolution operator 
U(t ) = exp(— itTi), namely 

W')) = 

In a coherent-state representation the evolution is effected by an integral kernel 


in the form 


K,(p", q", t*j p', q \ (') = (/, ,) 

q", i") = J K v (p", q", t"\ p', q’, t')$r){p', q r , t')dp'dq /2n . 


Clearly, K n depends strongly on the fiducial vector as does y n . 

Our goal in this paper is to formulate a universal propagator K{p" ,q" ,q' ,f'), a single 
function independent of any particular fiducial vector, which, nevertheless, has the property that 


Mp",q",t") = J K ( P ", 


q\ p', q', ,(p', q, t')dp dq' /2ir 


(1) 


holds just as before for any choice of fiducial vector. 

The functions A,, and K are qualitatively different as is clear from their behavior as t" —* t' . 
In particular 


Km = (p",q";r) \p' , q'-q). 


(2) 


which clearly retains a strong dependence on the fiducial vector. On the other hand, if (1) is to 
hold for any 77, we must require that 


Jm, K(p\ q", t"; p q\ t') = 2ir8(p" - p')6{q" - q') . 


( 3 ) 


Next let us turn our attention to a suitable differential equation satisfied by I< v and K. It is 
straightforward to see that 

1 ^) = \ P ^)' 

( 9 + *^p) b'T’lW = (P'WlQW’) 

hold quite independently of \q). Thus if H = H{P,Q) denotes the Hamiltonian it follows that 
Schrodinger’s equation takes the form 

Q 

= { p,iw\np,Q)m » 

d d 
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valid for any |p). 3 The propagators are also solutions of Schrodinger’s equation so it follows that 

where K # denotes either K n or K. What distinguishes which function is under consideration is 
the initial condition (at t = t 1 ) of the solution, namely, either (2) or (3). 

When K„ is under consideration, the operators -ij^ and q + i^ refer to a single degree of 

freedom made irreducible by confining attention to the subspace of L (R ,dpdq/ 27 t) spanned by 
ip v (p,q) for a fixed |r?) and for all |V>) € H. This restriction is implicit in K n because as t" -> f 
the resultant integral kernel (p", <?"; p|p', q'\ q) is a projection operator onto the subspace of an 
irreducible representation. 


1.2 The Universal Propagator 

In contrast to the former case, when the universal propagator K is under consideration the resul- 
tant Schrodinger equation (4) is interpreted as one appropriate to two degrees of freedom. In this 
view yi = q and y 2 = P denote two “coordinates”, and one is looking at the irreducible Schrodinger 
representation of a special class of two- variable Hamiltonians, ones where the classical Hamiltonian 
is restricted to have the form H c (p\,yi — P2), rather than the most general form H c (p\, p2,yi,y2)- 
In the case of K , and based on the interpretation described above, a standard phase-space path 
integral solution may be given for the universal propagator. In particular, and for a sufficiently 
wide class of Hamiltonians, it follows that 

Kip", q", t p\ q\ t') = M J j f [x ’ >+k< '-' H{k ' q - x)]dt VpVqT>kVx . 

Note that “x” and “A;” are “momenta” conjugate to the “coordinates” “p” and “7”, and also that 
the special form of the Hamiltonian has been used. In the standard phase-space path integral there 
is always one more ( k , x) pair of integrals compared to the (p, q) family, and the (k, x) integrals are 
unrestricted. This situation is made explicit in the regularized prescription for the path integral 
given, in standard notation, by 

I<{p",q",t"-,p', q\t') 

— lim 2 tt J J 

L-*oo 

L L 

x JJ dpidqt dk (+ idx e+ L / (27 t ) 2 , 
t = 1 (=0 

where p L +i,qL+i = p" , ?",Po, <7o = p',q', and where (L + l)e = ( t " - t') is held fixed. Let us first 
change the variables x t+ i — ► x /+ i + {qt+i + <7f)/2, followed by a second change x (+ l — > —x (+ i. 
The resultant regularized path integral reads 

K{p",q"X-,p',q',t') 
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_ 27T J • • • J e ^'=° [ 2 (?/+1+9*)(p/+1 -Pt)+k n £ (</<+! -9/)-x, + J. (p<+ 1 -p,)-eW(fc, + £ ,* <+ j. ) 


x JJ dp e dq ( JJ dk (+ idx t+ i. / (27r) 2 , 


<=i 


(=0 


or the formal path integral is given by 


I<{p",q",t"\p',q',t') = M J e i S [q *+ k i- xi, - H{ - k ' x ^ dt 'DpVqVkVx , (5) 


which is our final expression for the universal propagator in the present case. From this formula 
it is clear that the dependence on p" and p' is always of the form p" — p', and the dependence on 
q" and q' is always of the form q" — q' save for the universal phase factor \{q" + q'){p" — p'). In 
other words, 


K(p", q", t"; p\ 0 = F(p" - p\ - q\ (" - t‘)eW+W-r'-> 


for some function F. Of course, if 7i depends explicitly on time then F is not simply a function 
of the time difference t" — t' . 


2 Examples of the Universal Propagator 

2.1 Vanishing Hamiltonian 

Let us evaluate the universal propagator in three soluble examples. The simplest case is that of a 
vanishing Hamiltonian which leads to 

K(p",q",t"-,p',q',t') = M j jShv+^-^VpVqVkVx 

= M J Ff^SifiSifiVpVq 

= 2 n6(p"-p')6(q"-q'), 


where the normalization follows from the initial condition. Evidently this is the correct result. 


2.2 Free Particle 

The next case is the free particle where 7i(k,x) = k 2 /2m. In this event 

K(p", q", f p', q\ t') = M J e* / 1^'VpVqVkVx 

= jV J e i f {g * +m, ’ 2 M dt 6{p}VpVq 


j ^-Ap"-p'y 
i(t" - t ') ^ p } 


22 



2.3 Harmonic Oscillator 

The last case we consider is the harmonic oscillator where 7 = ( k 2 + u> 2 x 2 )/ 2. Now 


K(p",q"X;p',q',n 


= M J e i f w+kj ~ x *- {k * + “ ix 2 )/ 2 ]dt VpVq'Dk'Dx 
= M J e i f [q t + ( i 2 + *‘ 3 /u ' 1 )/ 2 ]di 'Dp'Dq 

= (2i)~ l csc(u>T/2) exp (i j ^(q" + q){p" -p') + J cot(wT/2) j^(p" “ p’Y + w (<?" 



( 6 ) 

where T = t" — V . Observe that this result is rather different from conventional propagators for 
harmonic oscillator Hamiltonians. Indeed, (6) is more like the propagator for a two-dimensional 
free particle in a uniform magnetic field. 4 This result also applies even when u — ► iw, or with a 
suitable limit, even when u> — ► 0 leading to the free particle solution. 

3 Propagation with the Universal Propagator 

In order to check our results for the universal propagator let us put them to the test. For ease of 
computation we choose as the initial state for our propagation the coherent-state overlap function 

{p'Y;p|p>?;»?) , 

and additionally we choose the fiducial vector |r/) to be the ground state of an oscillator with 
frequency f) for which (q\Q\q) — 0 = (p|P|p). In that case the initial state reads 

^{p',q’,t') = (p',g';p|p,9;p) 

= ex p(* {^(P +P)(?' - ~ ^[ fi-1 (p' ~P) 2 + W~?) 2 ]}) • 


3.1 Free Particle 

For the free particle case we need to compute ( T = t" — t') 

J exp { l -mT~\q" - q') 2 + l -(p" + p){q' - q) - ^[H -1 (p" - p) 2 + tt(q' ~ $) 2 ] j dq ' , 
which is readily found to be [p = (p" + p)/2, q = ( q " + q)/2,p* = p" — p,q* = q" — <7] 

Up", ?",<") 


\/rn 


+ ifIT/2 


exp 


{ 


im (Q, 2 Tq* 2 /i — Tp 2 + 2mpq*) p* 2 Cl(mq* — Tp) 2 | 


(m 2 + H 2 T 2 / 4) 


4H 4(m 2 + fl 2 T 2 /4) J 


( 7 ) 
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This result agrees with one obtained elsewhere 5 thereby establishing its validity. 

In addition, as fl — ► oo or fl — ♦ 0, Eq. (7), apart from a suitable scale factor, yields the result 
in the sharp q or sharp p representation, respectively. 3,5 In particular, consider 




-i 


2iriT 


which is the appropriate Schrodinger representation solution for the free particle. Likewise, con- 
sider 


= Iim 




mu / — c 

n-o 2\Arn 

= 8{p" — p)e~ 1 '*™ , 
which is the proper answer in momentum space for the free particle. 


3.2 Harmonic Oscillator 

Finally, let us consider the time evolution for the harmonic oscillator as given by 
(4 7 ri)- 1 csc(u;r/2) J exp (i fyq" 4 q'){p" - p’) + i C ot(u;T/2) \^{p" - p'f 4 u{q" - 9 ') 2 ] }) 

X exp j^(p' 4 p){q' -q)- jffr^p' - p) 2 + fl( 9 ' - 9 ) 2 ]} dp'dq ' , 

which is evaluated as [s = sin(u>r/2), c = cos(a;T'/2)]. 

u; 2 c 2 \ 

w 2 c 2 4 fl 2 s 2 ) 


0*J (p",q",t") = C(T) 2 exp iqp 


2 2 
U> Z S* 


J 7 2 C 2 + UJ 2 S 2 


+ ipq* 


x exp 
x exp 


VjOSC 


( 4 p 2 - ny 2 ) - 


iu>fl 2 sc 


where 


4(a> 2 c 2 4 fl 2 s 2 ) v ~ r " ■* 7 4(fl 2 c 2 4 u? 2 s 2 ) 

fl(cp* 4- 2ojsq) 2 flu > 2 (cq* — 2u r l sp) 2 

' 4(fl 2 c 2 4 u> 2 s 2 ) 4(u> 2 c 2 -f fl 2 s 2 ) 

C(T ) = cos ojT + i- ("—4 — ■) sin u>T . 

2 \Sl u J 


(49 2 -fl - 2 p* 2 ) 


(8) 


As in the free particle case, we can obtain the sharp q propagator by the same kind of limit, 
namely 

w.o = jim 

n^oo l V 7T 




2 7r z sinter 


exp { l cos “ T(, " ! + - 2 ’"«] } ■ 
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which, of course, is the standard result. Likewise, the sharp p propagator is given by 


1 


\/ 27 ria; sin uT ^ 1 2u> sin u>T 


{■ 


[cos u>T(p" 2 + p 2 ) - 2p"p 


which, again, is the standard result. , ., , , , 0 _ 

We may also observe that the harmonic oscillator evolution simplifies considerably when il - w 

In that case 1 


V>„(p", 9 ", *") = exp —(9 w p w - qp ) + - pV) - ^(p" - Pt) 2 4 (?" 9*0 » 


where 

= gcos(uiT) -f uT^psin^T) , 
pj = p cos(air) - u><7 sin(u>T) , 

evolution equations that are seen to follow the classical solution. 

3.3 Generalization 

Although we have only shown that a limited set of fiducial vectors are correctly propagated by the 
universal propagator, it should be fairly clear that the stated properties of the universal propagator 
hold true. Indeed, the general case may be discussed by considering as initial condition 

(p, q; rj\e- ip ' Q e iq ' P \p, q; rj) = e~ <fr V + p,q' + 95 q\P, 9i v) 

= e '(P4'-?P')+«P ( ?- | I ) ^p' i q\p-p,q- q) 

for just a single |p>, say a Gaussian with 0 = 1. Then a suitable superposition over p,q leads 
to any fiducial vector of interest, while a second and independent suitable superposition over p, q 
leads to any initial state |^>) of interest. 


4 Classical Limit 

Although the universal propagator has been derived by identifying the relevant Schrodinger equa- 
tion as one for two degrees of freedom, it should nevertheless be true that the classical limit refers 
to a single degree of freedom. This is possible, in the present case, because of the limited form ot 
the quantum or classical Hamiltonian. 

Recall, under standard assumptions, that the classical action for a conventional coherent state 
path integral reads, in the limit h — ► 0, as 

/ = J\pq- ft(p> q)] dt • 
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Extremal variation of this expression holding the end points fixed leads to the usual Hamiltonian 
equations of motion, 


q = d'H(p,q)/dp, 

P = ~9'H(p,q)/dq, 

appropriate to a single degree of freedom. Let us denote a generic solution of these equations bv 
<? c (f) and p c (t). 

Before proceeding it is important to reexamine the “standard assumptions” that lead to this 
result. For finite h the expression that represents the classical Hamiltonian in a coherent state 
path integral is traditionally given by 

H{p,q) = (p,qiv\'H(P,Q)\p,q-,Ti) 

= (q\n(P + p,Q + q)\q) . 

Normally, one restricts fo) so that (q\ Q\q) = 0 = (r}\P\rj), and ( q\Q 2 \q ) -+ 0 and (t/JP 2 ^) ^ 0 as 
h — ► 0. In this case 

limtf(p, 9 ) = H{p,q). 

However, in the present paper we want to deal with more general fiducial vectors \q) such that 

( vlQlv ) = qr,, 

(t)\P\t)) = p„, 

are arbitrary real variables. We still insist on vanishing dispersion as h — ► 0, namely, that 

(v\(Q - qrflv) -o, 

(v\(P - Pv) 2 \q) -o, 

as h -» 0. This more general situation leads to the result 

lim H (p, q) = H(p + q + 9 „) 

as the representative of the classical Hamiltonian. 

In this more general case the classical action appropriate to the coherent state path integral 
becomes & 

I ~ J [(p "b Pv)q — q^p ~ %(p + Pi?» q + ?»))]<& • 

In this expression p = p(t) and q = q(t), while p^ and q v are time-independent constants. The term 
f(Pf?q - qvP) dt = Pv(q" - q f ) - qv(p" - p') is a pure surface term and will not affect the equations 
of motion; it could be eliminated simply by a phase change of the coherent states. Extremal 
variation leads to the equations of motion 

q = dHip + Prj^ + q^/dp 
P = -dft(p + Pr,,q + qv)/dq, 
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which have as their solutions 


q(t) — <]c(t) <1 t) > 

p(t) = p c (t) - p„ , 

where q c (t ) and p c (t) denote a generic solution of Hamilton’s equations when q n = p v = 0, as 
discussed above. 

Finally, we note that although the dispersion of \tj) vanishes as h -» 0, the generally nonva- 
nishing values of q„ and p„ are vestiges of the coherent-state representation induced by |r/) that 
remain even after h — ► 0. 

4.1 Classical Limit of the Universal Propagator 

In the case of the universal propagator the expression that serves as the classical action is identified 
as [cf. (5)] 

I — J [qp + kq — xp — 'H(k,x)]dt . (9) 

Extremal variation of this expression holding the end points fixed leads to the set of equations 

q = X, 

p = ic, 

q = dH(k,x)/dk , 
p = —d'H(k,x)/dx . 

Consequently 

x = d'H(k,x)/dk , 
k = —d7i(k,x)/dx , 

which show that (k,x) satisfy exactly the same equations of motion as do ( p c ,q c ) in the usual 
classical theory. Thus we may identify the solution k(t),x(t) with p c (t),qc(t)- In addition, we 
have 

q(t) = 9e(0 “ c? * 

p(0 = Pc{t) ~ C p i 

where c, and Cp denote two arbitrary integration constants. Among all possible values of c q and 
Cp are those that coincide with q n and p n for a general \rj). 

Thus we find that the set of solutions of the universal classical equations of motion appropriate 
to the universal propagator includes every possible solution of the classical equations of motion 
appropriate to the most general coherent-state propagator (with 1^) having vanishing dispersion 
as h -+ 0). Not only does the quantum dynamics (universal propagator) correctly evolve the state 
vectors in a canonical coherent-state representation for a general |t/), but the classical dynamics 
(universal classical equations of motion) correctly evolves the classical phase space points accord- 
ing to the coherent-state induced classical equations imprinted with arbitrary values of the only 
remnant of the fiducial vector after h -> 0, namely its average coordinate and momentum values. 
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5 Extension to Other Coherent States 


We observe that the procedure to introduce a Schrodinger equation and a path integral solution 
for the universal propagator applies for other sets of coherent states, such as the spin coherent 
states, the affine [or SU(1,1)] coherent states, etc. 6 In each of these cases it becomes possible to 
introduce an appropriate universal propagator just by following the procedure we have given for 
the canonical coherent state case. 
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Abstract 

We use this simple example to show how the formalism of Moyal works when it is applied 
to systems of identical particles. The symmetric and antisymmetric Moyal propagators are 
evaluated for this case; from them, the correct energy levels of energy are obtained, as well 
as the Wigner functions for the symmetric and antisymmetric states of the two identical 
particle system. Finally, the solution of the Bloch equation is straightforwardly obtained 
from the expressions of the Moyal propagators. 


1 Phase-space Q M formalism 

The original ideas of this approach to Q M are due to Weyl [1], Wigner [2] and Moyal [3]. States 
and observables are no longer operators on a Hilbert space but functions on an adequate phase 
space. The Weyl mapping relates both formalisms: given a function / defined over the phase 
space i? 2n , the corresponding operator F is given by 


p = /„„/(“) n < u )^ u = (s-P)- 


( 2ith) 

Reciprocally, given an operator A the associated function in the phase space is 

Mu) = tr {in(u)} = W-' (A). 

As we can see, a central role is played by the “Grossman-Royer” operators [4, 5]: 

^( 2 q - rj). 


[ff(q, p)'P](77) = 2 n exp 


LX p(t? ~ q) 


( 1 ) 


(2) 


( 3 ) 


'On leave of absence, Departamento de Fisica Teorica, Universidad de Valladolid, Spain. 
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Cs ti i 


The twisted product of two functions is defined as the non-commutative operation that corre- 
sponds to the product of operators: 


(/x 9 )(u) = W-'(WU) W{gj) 

= (Spt /(v|,w ' ,p 

where the matrix J is simply 

0 1 

l -I 0 


2 i 

— (uJv + vjw + wJu) 
n 


dv dw, 


(4) 


(5) 


being I the n-dimensional identity matrix. 

In the Schrodinger representation of quantum mechanics, the information about dynamics 
is contained in the evolution operator U(t). Its counterpart in this formalism is the “Moyal 
propagator”, defined as 

Z(u,t) = W~ l (u(t)). (6) 

It verifies Schrodinger equation: 

ih w = HxE - < 7 > 

The Fourier transform of this function with respect to t gives the spectral projections parametrized 
by E : 

nu ' E) = ^fil R E(n ' t)e ' E, "' dt (8) 

If the Hamiltonian is time independent, the support on E of T coincides with the spectrum of 
H [6]. If E 0 belongs to the discrete spectrum of H, T(u, E 0 ) is, but for a constant factor, the 
Wigner function of the orthogonal projector into the proper subspace Eq [6]: 


Wv,(q,p) = (MM) = cdidr Jr- e ' py/ V(q + y/2)^(q - y/2) dy. 


( 27 vh) 


(9) 


2 Phase-space Q M formalism for identical particles 

In the standard formalism of quantum mechanics, to deal with a system of N identical particles, 
we introduce a superselection rule: the space of physical states is a closed subspace of the initial 
Hilbert space. The Hilbert space is splitted [7] 

H = ( 10 ) 

where ?{+ is the Hilbert space of the wave functions symmetric under the exchange of any two 
particles and the Hilbert space of the antisymmetric functions. The functions in 7i w have no 
symmetry of this kind. The orthogonal projectors are given by 

r+ = -& Z P’ • c- = 4t (ii) 

tfCPjV ’ <7 €Pat 
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where 


( 12 ) 


(■^ > uV , )(^'l i • ■ ■ ) ^n) — <7(1)1 • ■ • 1 -^<j(A r ))- 

for any <7 in the group P,v of permutations of N elements; 7r(<r) is the parity of a. 

If B is either an observable or an state in 7 ~t, the corresponding operators for a system of A 
fermions or bosons are: 

P_BP_ and P+BP+. (13) 

If B is invariant under the exchange of particles, we have 

P ± BP ± = BP ± = P±B. (14) 

We use the Weyl transformation to translate these ideas into the language of phase space. 
Therefore, the function for an state or observable when we consider N bosons or fermions is 

W~\P ± BP ± ) = W~\P ± ) x W~\B) x W~\P ± ). (15) 

Due to the fact that the Weyl map is linear, all we need is the function for any permutation 
a. As a can be written as the product of cyclic permutations with no common elements [8], it 
is enough to compute the function corresponding to such a cycle. If we consider a general cycle 
<7 = (1,2,3, ...,M) we get: 

a = 2 (M - 1)n exp{-^ Y {-l) k+l u*Juil, M odd; (16) 

{ h k=i 1>k ) 

(o i M 1 

a = (2 m_1 7rft) n (5(u 1 - u 2 + . . • - u w ) exp < -y E (-l) fc+ ' u,J Ul , M even. (17) 

[ n k=V,l>k J 

As an example, for a two cycle that exchanges the particles i and j we have: 

&ij{ Ui,...,u/v) = = (27rfi) n 6(u, -Uj), (18) 

and it can be checked that 

{crij x p x a t] ) (. . . u„ . . . Up . . .) = p{. Uj, . u,-, . . .). (19) 

The functions corresponding to the orthogonal projectors for a system of two onedimensional 
particles are 


p+(U!,U 2 ) = 

hi + 27r/i<5(ui - u 2 )), 

(20) 

P— ( Ui,U 2 ) = 

1(1 - 27rfi6(u! - u 2 )). 

(21) 
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3 Two onedimensional identical particles under an ex- 
ternal oscillatory potential 

Along the present section, we intend to presenting an example of particular interest in order to 
illustrate the preceeding discussion. We shall study the behavior of a two onedimensional particle 
system subjected to oscillatory forces of the same frequency. If we do not take into account the 
identity of the particles, the Hamiltonian will be simply: 


H( u„u 2 ) = H( u.) + H( u 2 ) = + p\) + + ,|). 

The corresponding Moyal propagator has been already evaluated [6], and is: 

S(ut,u 2 ) = H(u,)H(u 2 ) = g exp j-^(ff( u i) + W(u 2 ))tan yj. 


( 22 ) 


(23) 


Let us now introduce the statistics. As 7 /(ui,u 2 ) is invariant under permutations of the tw r o 
particles, the Hamiltonian for our system of two identical particles is: 


H ± (u l ,u 2 ) — (Hx p± )(uj , u 2 ) = ~ (ui ) + H( u 2 ) ^ 27rft 3 | 


' ^(u! — u 2 ) mu> 2 S( u x — u 2 )^ 

2m(q 1 - q 2 ) 2 + 2 (p, - p 2 ) 2 j 


We see that, after symmetrization or antisymmetrization, the Hamiltonian on phase space of 
our system is not longer (22) but (24). Equation (24) includes (22) plus an extra term. From this 
term results an extra potential, due to the introduction of the statistics, which has a quite different 
action depending whether the particles are fermions or bosons. In the first case, this potential is 
preceeded by a plus sign and, therefore, it is equivalent to a delta barrier preventing that qi = q 2 
and pi = p 2 . This already suggests that both particles cannot remain in the same state and, 
hence, that they fulfill the Pauli principle. This idea will be confirmed by our calculations for the 
lowest energy levels. On the contrary, if the particles are bosons the extra term has a minus sign 
and, consequently, it represents the apparition of a delta well. This delta well would rather favor 
the presence of particles in the same quantum state. In a clear opposition to the case of fermions, 
no exclusion principle can exist here. 

The symmetrized Moyal propagator is obtained in a similar way: 


-±(ui,u 2 ) 


= - 1 

2 \ cos 2 f 


f-2i(H(u l ) + H(u 2 )) i ut' 

exp i — r 1 — — tan — 

F V hu 2 


±{-i) 


it exp 1 

+ U 2) tan f] 

f ex Pi 

[^(u,-u 2 )cotf] 

LI 

l ) 

cos y 


sinf ( 


(25) 


Comparing with (23) we see that there is also an extra term due to the statistics. 
The spectral projections are obtained from (8) and (25); in this case we obtain 


r±(u 1 ,u 2 , J E) = 2e- 2 ^ H{Ul)+H ^ /hw f^{-l) k 6(E- hu(k + 1)) 


k 

x £ 

n=0 


2/f( Ul + u 2 )\ (2H{ui-n 2 f 
hu ) k ~ n { hu j 


.(26) 
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From here, the well know energy levels are obtained for the fermionic and bosonic cases. Let us 
notice the coefficient of 6(E — hu), that vanishes for fermions but not for bosons. 

We can evaluate the Wigner functions corresponding to states of two particles, in states i and 
j. Let us write those functions as W,“, the corresponding to the antisymmetric state, and W,* the 
associated to the symmetric state. We then have: 

r + (u„u 2 ,£) = (2rt) 2 [ W£(q,p)<(£-*w) + Wi(q,p)*(£-2Aw) 

+ (W,',(q,p) + W* 2 (q,p)) 6 (E - 3 M + •••], (27) 

r_(u,,u 2 ,£) = (2>rft) 2 ( w;,(q,p)«(£- 2*u.) + W 0 “ 2 (q,p)«(£ - 3M +••■]. (28) 

The coefficients of the 8 are the Wigner functions of the orthogonal projector on the corresponding 
eigenspaces. 


To finish, let us solve the Bloch equation, that is, let us find the Wigner function corresponding 
to the density matrix of the canonical ensemble for the system we are considering. Bloch equation 
reads simply: 

^ = -H x O = -fi x //, (S=l/kT, (29) 

that is, it is Schrodinger’s equation with the change t —* —ihft. 

But, as we already know the form of the Moyal propagator, we can write inmediately the 
solution for Q(u 1 ,u 2 ,^) by making the change t —* —ih/3 in E ± (u 1 ,u 2 ,<). We get: 


ft±(«i,u 2 ,0) 


= if 


2 (rosh 2 
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ex pfe^( u i + u 2 )tanh^} exp{^ff(u! - u 2 )coth^} 


cosh 


sinh^ 


(30) 

.(31) 


After integration of f)(uj,u 2 ,/?) over the phase space, we get 

exp(±fiw/?/2) 


Z ± ((3) = 


8 cosh(/tu;/?/2) sinh 2 (/iu;/?/2) 


(32) 


From this partition function, we can obtain the thermodynamical quantities, for example the 
internal energy, the free energy and the entropy 


E±m = 
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2 2 J 
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+ log 8 } T 
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(33) 

(34) 

— In sl . (35) 


Notice that the entropy is the same in both cases (bosonic and fermionic). 
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